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Abstract 

This paper offers a model for incoherent scatter signal spectra with- 
out averaging the received signal over sounding runs (realizations). The 
model is based on the existent theory of radio waves single scattering from 
the medium dielectric permittivity irregularities, and on the existing ki- 
netic theory of the plasma thermal irregularities. The proposed model 
is obtained for the case of monostatic sounding. The model shows that 
the main contribution to the received signal is made by ion-acoustic waves 
caused by certain spatial harmonics of the ions number phase density. The 
model explains the width and form of the signal spectrum by macroscopic 
characteristics of the medium, and its fine 'peaked' structure by character- 
istics of the ions number phase density. The notion of the weight volume 
is introduced to define the domain of wave vectors and velocities in the 
spatial spectrum of the ions number phase density which makes the main 
contribution to the formation of the scattered signal. This weight volume 
depends on the antenna pattern, the form of the sounding signal, and on 
the time window of spectral processing, as well as on ionospheric plasma 
macroscopic parameters: electron and ion temperatures, ion composition, 
and the drift velocity. Within the context of additional assumption about 
the ions number phase density, the proposed model was tested by the data 
from the Irkutsk incoherent scatter radar. The test showed a good fit of 
the model to experiment. 



1 Introduction 

One of the remote probing techniques for the ionosphere is the method of radio 
waves incoherent scatter. The method is based on the scattering of radio waves 
from ionospheric plasma dielectric permittivity irregularities [Evans, 1969]. Fur- 
thermore, two different experimental configurations are involved: monostatic 
(where the receive and transmit antennas are combined) and bistatic (where 
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these antennas arc spaced). In actual practice, it is customary to use the mono- 
static configuration. Ionospheric plasma parameters (ion composition, drift ve- 
locity, electron and ion temperatures, and electron density) in this case are 
determined from the scattered signal received after completion of the radiated 
pulse. The spectral power of the received signal, averaged over sounding runs 
('realizations'), is related (assuming that such an averaging is equivalent to 
statistical averaging) to the mean spectral density of dielectric permittivity ir- 
regularities by the radar equation [Tatarsky, 1969]. The connection of the dielec- 
tric permittivity irregularities spectral density with mean statistical parameters 
of the medium is usually determined in terms of kinetic theory [Clemow and 
Dougherty, lQm\ Sheffield, 1975; Kofman,imi]. 

The location and size of the ionospheric region that makes a contribution 
to the scattered signal (sounding volume) is determined by the antenna beam 
shape, the sounding radio pulse, and by the time window of spectral processing 
[Suni et al., 1989]. The shape of the sounding volume determines also the 
method's spectral resolution, the accuracy to which the mean spectral density 
of dielectric permittivity is determined (which, in turn, affects the determination 
accuracy of macroscopic ionospheric parameters: electron and ion temperatures, 
the drift velocity, and electron density). The number of realizations, over which 
the received signal spectral power is averaged, determines the method's time 
resolution, i.e. its ability to keep track (based on measurements) of fast changes 
of macroscopic parameters of ionospheric plasma. 

Currently most incoherent scatter radars have accumulated extensive sets of 
the scattered signal individual realizations (private communications of P.Erickson 
(Millstone Hill), V.Lysenko (Kharkov IS radar), and G.Wannberg (EISCAT)). 
Therefore, attempts are made to analyze the realizations from different methods 
which differ from a standard averaging by their sounding runs. Basically, these 
methods imply looking for small scatterers making the main contribution to 
the scattered signal. This method is good for analyzing signals scattered from 
meteors and their traces [Pellinen-Wannberg, 1998]; however, it is insufficiently 
substantiated for describing the scattering in the ionosphere. 

In the work there were used the experimental data obtained with Irkutsk 
Incoherent Scatter radar. The radar is located at 52° A/', W4PE, it has sounding 
frequency 152-160 MHz and peak power 3MW. High signal-to-noice ratio during 
the experiments under investigation ( S/N ^ 10 ) allows us to neglect the noice 
effects when analyzing the signal received. 

The technique of the incoherent scatter signal processing in Irkutsk IS radar 
is the following. For each single realization of received signal we calculate spec- 
trum in time window with width equal to the sounding signal duration and with 
delay corresponding to the radar range to the sounding volume investigated. 
The sounding signal we use in this experiment is a radiopulse with duration 
800 mks. The repeating frequency approximately 25 Hz. Averaging over the 
1000 realizations corresponds to 3% dispersium of the overaged spectrum rela- 
tive to its mathematical expectation. The reason of using such a simple pulse 
is to investigate the fine structure of the single (unaveraged) spectrums in this 
simpliest case. 
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Figure 1: The spectral power of three different reahzations (A-C), averaged 
over 1000 reahzations spectral power of the scattered signal (D) and spectral 
power of the sounder signal envelope (E) as deduced using the data from the 
Irkutsk incoherent scatter radar. 

Figure |^ exemplifies the mean spectral power of the scattered signal and 
its separate realizations, based on the data from the Irkutsk Incoherent Scatter 
radar. It is evident from the figure that the spectral power of the scattered 
signal in an individual realization (Figure |^(b-d)) differs drastically from that 
averaged over realizations (Figure 0(a)); therefore, existing model of the incoher- 
ently scattered signal, based on averaging over sounding runs, are inapplicable 
for its interpretation. For that reason, development of new models of the scat- 
tered signal for analyzing its separate realizations without averaging them is 
important from the theoretical and practical standpoint. 

Sometimes it is useful to suppose that incoherent scattering signal is a ran- 
dom gaussian one [for example, Farley, 1969; Zhou, 1999]. But, it is well known 
that the signal received is a detirministic function of ionospheric dielectric per- 
mittivity e and is fully determined in first approximation by the Bern's formula 
(in one or another its form [Ishimaru, 1978; Berngardt and Potekhin, 2000]), this 
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relation could be called as a radar equation for signals [Berngardt and Potekhin, 
2000]. 

The dielectric permittivity irregularities also could be supposed as a random 
functions, but they are deterministic functional of some other functions (in 
case of uncollisional unmagnetized palsma with one ions type those functions 
are phase density of the ions and electrons as functions of velocity, location and 
time, ion composition and temperatures of the ions and electrons, this functional 
dependence is determined by the Landau's solution [Landau, 1946]). 

If one could determine all these unknown functions, the received signal shape 
in single realization will be fully determined, and could be analyzed without us- 
ing any statistical methods. Such an approach, for example, is used in radioa- 
coustical technique of the athmosphere sounding when the delectric permittivity 
irregularities (by which the radisignal is scattered) are generated by the acous- 
tical wave [Kalistratova and Kon, 1985]. 

The statistical proporties of the single realizations are showed at Figure ||. 
From this figure it becomes clear that the unaveraged spectrum has the fine 
structure - it consists from a number of peaks with approximately 1.5KHz width 
(and this width very slightly depends on frequency), which could be character- 
ized by the peak amplitude (amplitude at the maximum of the peak) and peak 
appearence (number of realizations in which there is a peak maximum at given 
frequency) at the given frequency, and those properties distributions are not 
gaussian ones but have double peaked structure and located in the same band 
with incoherent scattering average spectral power. 

This fact allows as to suppose that not only average spectral power of the 
received signal depends on ionosperical parameters, but the fine structure of 
non-averaged spectra too. 

At first, it is neccessary to understand quilitatively, what information one 
could obtain from one realization of the IS signal. It is well known, that af- 
ter any statistical processing of a function a part of the information is loosed 
irreversibly (for example, when one calucaltes the first n statistical moments, 
all the rest moments, starting with n-l-1 are still unknown). That is why, if 
the statistical characteristics of the realizations (mean spectral power or corre- 
lation function) are depend on the ion and electron temperatures and the ion 
composition then single realization must depend on all those parameters and on 
some new 'additional' parameters. It is clear that to determine temperatures 
and ion composition from averaged signal parameters is much easier than from 
single realization (because the second one includes additional parameters), and 
we can use the ones obtained from mean spectral power, with necessary spa- 
tial and spectral resolution, using different techniques, for example alternating 
codes [Lehtinen, 1986]. But the new 'additional' parameters can be determined 
from single realizations only. 

The aim of this paper is to find out the functional dependence of single 
realization spectrum on all the parameters, including well known (temperatures 
and ion composition) and new ones, which could describe the single realizations 
spectrum properties. For this propose we will use for analysis only signals with 
high signal to noise ratio (more than 10), because in this case the noice effects 
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could be neglected and the received signal could be supposed as only IS signal 
without presence any noice. 



2 Initial expressions 

To analyze the individual realizations of the scattered signal, it is necessary to 
have a convenient expression relating the spectrum of the scattered signal to the 
space-time spectrum of dielectric permittivity irregularities without averaging 
over realizations. Such an expression for a monostatic experimental configu- 
ration was obtained and analyzed in [Berngardt and Potekhin, 2000]. It holds 
true in the far zone of the receive-transmit antenna and, within constant factors 
(unimportant for a subsequent discussion), is 

u{uj) = [ H{uj -v,k-2ko- ~t)dvd~t. (1) 

Here t(y^ fc ) - is the space-time spectrum of dielectric permittivity irregu- 
larities; H{Lo,k) « / H{t,r)e-''^'^^+'''^'>drdt = / o{t)a{t - 2r / c)e-^'^'^*+^''Utdr / r 
- is the narrow-band weight function; a{t),o{t) - are, respectively, the sounder 
signal envelope and the time window of spectral processing; g{r) - is the antenna 
factor which is the product of the antenna patterns by reception and transmis- 
sion; r = V'/r - is a unit vector in a given direction; fco - the wave number of 
the sounding wave; c -is the velocity of light. 

Suppose that sounding signal and receiving window of spectral processing 
are located in time near the moments t = Ti — Tq and t — Ti respective! 
and theirs carriers do not intersects (this is the one of the radar equation (| 
obtaining conditions [Berngardt and Potekhin, 2000]). In this case the carrier 
of the weight function H{t, r) is located near the t = Ti; r = Tqc/2. By going in 
equation (|lj) to the spectrums calcualated relative to the weight volume center 
(to remove the oscillated multipliers under integral), we obtain (neglecting to 
the unessentional multiplier): 

u{uj)^ jH^{y^~v,k-2ka-fJ^_^. (2) 
e{v, ~t- Ti - To/2, -kToc/2)dvdli 

where Hi{uj, k) H{uj, fc)e*'='^o^/2 _ low oscillating part of the H{uj, k), cor- 
responding to its calculation relative to the center of the weight volume H{t, r); 

and ?(i^. It; T, it) — e{i', li)e'^^'^'^^ ^ - is a time-spatial spectrum of dielectric 
permittivity irregularities calculated relative to the point t = T, = 7?. 

In accordance with [Shejfield,1975; Clemmow and Dougherty,1969; and Akhieszer 
et a/., 1974], assume that the spectrum of small-scale dielectric permittivity ir- 
regularities is determined by the Landau solution [Landau,19A9]. Then the 
low-frequency (ion-acoustic) part of the irregularities spectrum in a statistically 
homogeneous, unmagnetized, coUisionless ionospheric plasma with one sort of 
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ions is determined by plasma macroscopic parameters (electron and ion tem- 
peratures, ion composition, and drift velocity), and by unknown conditions in 
the moment T related to which this spectrum is calculated - the ions number 
phase density in a six-dimensional phase space of velocities and positions of 
particles. It is known that the dielectric permittivity irregularities spectrum at 
large wave numbers of the sounding wave fco > Wat/c (where wjv is plasma fre- 
quency) , is proportional to the electron density irregularities spectrum [Landau 
and Lifshitz,1982, par. 78]: 
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47ro 

e{u,,k;T) = --^n,,{u,k;T), 

which is given by the expression (for example, [Sheffield, 1975, sect. 6]): 

-j^ Ge(w,"fc*) f exp{itT')fail^,lf]T) 

nei{uj, k;T) = =;— / =f -drdv, (3) 

eii{u!, k) J ui — klf — ij 

where e\\{uj,lc) = (l + Ge{ijJ,~^) + Gi{u},~^)^ - is longitudinal dielectric 
permittivity; wave number k should be small enought to wave length be smaller 
than Debye length (Solpiter approximation). Most part of IS radars have the 
sounding frequencies(50-1000 MHz) within these limitations. 



Ge,i(a;, /c) = -2 J =^ dlf; (4) 

™e,i/C UJ — k V — I'y 

fe,io{'v^), ne,io- are equilibrium distribution functions of the electrons and 
ions velocity and their densities; me,i, qe,i- are the mass and charges of electrons 
and ions, respectively; 

N 

fa{r,lt;T)= E ^{'r'-n{Tm-if-v^{T))-U{ir) (5) 

i = i 

- the ions number phase density in a six-dimensional phase space of velocities 
and positions of particles (the ions number phase density, INPD, at t ~ T), 
with the summation made over all ions. Generally equilibriTim distribution 
functions fe,io{~^) are taken to be Maxwellian, with the temperatures and 
Tj for electrons and ions, respectively, and in the absence of a drift they are 

/e,iO = exp{-{v/VT^^^f)/(^Vlf'^ 
1/2 

where vt^ . = (2fcT'e^i/me_i) ' stands for the thermal velocities of electrons and 
ions, respectively. Then the functions Ge,i{(jJ, k ) have the well-known analytical 
expression, for example [Sheffield,1975]: 
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(6) 



The physical meaning of the expression (|3|) is as foUows: the position and ve- 
locity of each ion at the moment T are determined by the INPD fii{~r^ , ~v^', T), 
and the dielectric permittivity irregularities l^{uj, k ; T) are determined by ion- 
acoustic oscillations of plasma under the action of such initial conditions. 



3 Traditional processing of the incoherent scat- 
ter signal, and characteristics of its separate 
realizations 

Traditionally, the incoherent scattered signal is processed in the following way. 
A set of the scattered signal spectra (|l|) is used to obtain its spectral power 
averaged over realizations. By assuming that an averaging over the realizations 
is equivalent to a statistical averaging, and also by assuming a Maxwellian 
distribution of the INPD fni'r' , T), one can obtain the following expression 
for the mean spectral power of the scattered signal [Suni et a/., 1989]: 



const 2nqi 
i?2 1^ 



F{uj - v) 



Ge{v,2ko) 



2ko 



(7) 



where F(lo) is the 'smearing' function determined by the spectrum of the 
sounder signal and the spectral processing time window; and <> is averaging 
over realizations. 

The frequency dependence of the scattered signal mean spectral power 
under usual ionospheric conditions has a typical 'two-hump' form (Figure |l|(a)) 
[Evans,1969]. From the scattered signal mean spectral power (Q) it is possible 
to determine the electron Te and ion Ti temperatures and the drift velocity 
ifo involved in a familiar way in the functions 2ko), Ge(w, 2fco) and fio(v) 
[Sheffield,m5]. 

In single realizations, however, the scattered signal spectral power differs 
essentially from the mean spectral power. Figure |l| presents the scattered sig- 
nalspectral power in three consecutive realizations (Figure ^(A-C)), the spectral 
power averaged over 1000 realizations (Figure |^(D)), and the spectral power of 
the sounder signal envelope (Figure |l|(E)). From Figure |l| it is evident that the 
non averaged spectral power of the incoherent scatter signal (Figure |l|(A-C)) has 
a typical 'peaked' form, the width of peaks is larger than that of the sounder 
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Figure 2: Statistical properties of single spectrum realizations (averaged over 
1000 realizations) as functions of frequency. Mean peak amplitude (A), mean 
peak width(B) and mean appearence of the peaks (C) for different frequencies 
of received signal. For comparison there is showed the IS spectral power for this 
data (D). 
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signal spectrum, and the peaks themselves are concentrated in the band of the 
mean signal spectral power. In the case of an averaging over realizations, such 
a peaked structure transforms to a typical smooth two-hump structure (Fig- 
ure|l|(D)). 

4 Model of single realizations of incoherent scat- 
ter signals 

4.1 Structure of one realization of the incoherent scatter 
signal spectrum 

To obtain a model of the incoherent scatter signal we substitute the Landau's 
expression for dielectric permittivity irregularities (|^) into the expression for the 
scattered signal spectrum (||). Using in the spatial spectrum k ,lt]T) 
of the ions number phase density /ii(T*, 1?; T), and upon interchanging the 
order of integration, we obtain the one-realization model for the scattered signal 
spectrum: 

uiu;)^ [ KiLj,t,v\\)Fait,v\\;Ti-To/2)dlidv\\, (8) 



= a^.k) d.. (9) 

Here 

F^^(^ ,v\\-T) = j Ui{l^,ir-T)5{tir -kv\\)dir (10) 

is unknown function we want to determine from experiment and has a form 
similar to the Radon transform of the function /ii( fc , it; T). A kinetic function 
(showed at Figure ||) 

f(i^,t)=Ge(i^,t)/e||(j.,fc*) (11) 

is determined by macroscopic parameters of ionospheric plasma Tg^j and vq; 
these parameters can be determined, for example, from measurements of the 
mean spectral power of the received signal (^ . 

Thus the kernel Kiuj, k ,v\\) is completely determined by the sounder signal, 
the receiving window, and by macroscopic characteristics of ionospheric plasma. 
The expression (^) clearly shows the meaning of the kernel K{uj, fc , Wjj): it deter- 
mines the selective properties of the model, i.e. the possibilities of determining 
the unknown function Fii{ k ^v\\]T) from the measured u{oj). Hence it can be 
termed the weight volume in the space ( , W||), or ambiguity function. Since the 
function Hi{uj, k) is a narrow-band one, with its carrier concentrated near zero, 
the function of indefiniteness K{uj, k,v\\) has also a limited carrier in k near 
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k — 2kQ. The possibilities of determining the unknown function Fii{ k ,U||;T) 

dependence on the wave vector directions k are determined by the product of 
the kinetic function ^(w, k ) and the antenna beam g{—k). 

According to the resulting model (||), the form of scattered signal single 
spectrum is determined both by a determinate component (the weight volume 
-fC(w, fc, W||)), and by a random (i.e. dependent on time by the unknown way) 
component. A random component is the function Fii{ k , T) determined by 
the spartial harmonics packet of the INPD fii{ k ,lf]T) with wave numbers 
k, concentrated near 2fco and calculated relative to the moment t = T. The 
moment T is determined by the moments of the sounding signal transmitting 
and spectral processing receiving window location, and corresponds to the mid- 
dle moment between them T = Ti — To/2. The weight volume K{uj,k^v\\) 
determines the parameters of this wave packet, the region of wave vectors and 
velocities in Fn (k ,v\\\T) which make the main contribution to the scattered 
signal at a particular frequency lo. 



4.2 Qualitative properties of the weight volume 

In this experiment the spectra sounder signal a(w) and the receiving window 
o{uj) are selected such as they are sufhciently narrow-band ones (IKHz), in com- 
parison with the functions e||(w, 2fco), Geito, 2fco) and fio{v), in order to improve 
the accuracy of their determination from experimental results. Therefore, the 
weight function i?i(w, fc) ~ o(u; — ck/2)a{ck/2) can also be considered narrow- 
band from both arguments as compared with the kinetic function ^(v, k ) from 
corresponding arguments. In this case we can approximate the weight volume 
K{uj, 1 as: 

£,{uj,2kok)g{-k) [ Hi{uj ~ v,k ~ 2ko ~ v / c) 
K{u;,k,v^^)^ p J (,_fc,„_,^) (12) 

The function Hi{Lu,k) is concentrated near lo — 0, k — [Berngardt and 
Potekhin,2QQQ\. The width of this function from arguments (oj, fc) is Aw = 
(Awo + Awo), Afc = (Awa + Au:o)/c, where Awq, AoJo, is the width of bands of 
sounder signal spectra and of the spectral processing window, respectively. In 
this experiment on ionospheric sounding by the incoherent scatter method, they 
have the order AcUa.o ~ W^sec^^. The function ^(i^, fc ) at a fixed k = 2ko = 
6.28to~^ for a typical ionospheric plasma with 0+ ions and = Ti ~ 1500K 
is presented in Figure ||. The figure shows that the function ^{ly, k ) does varies 
smoothly on the characteristic size of the weight function Hi(uj, k) carrier in u) 
which is in this case has the order of AH. ~ 10* sec~^ (corresponds to a sounding 
by the impulse radio signal of a duration of 1 millisecond) . 

Assuming that the envelope of the sounder pulse a{t) and the receiving win- 
dow o(t) have an identical Gaussian-like spectrum: a{u!) — o(w) — exp{~(uj/ Alu)'^) 
we obtain the function Hi(uj, k) of the form 

Hiiiv, k) = W2)/A^)^g-(W(2A^))^ (13) 
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Figure 3: Kinetic function S,{i',k) as a function of u at a, fixed k = 2ko = 
6.28m~^ for a typical ionospheric plasma with 0+ ions and Tg = Tj = ISOO-fT. 
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Upon substituting (|13|) into (|12|) and rather unwieldy calculations, similar to 
[Landau, 1946], we obtain the following expression for the scattered signal spec- 
trum (Bf): 



u{lu) = j K{uj,k,v\{)F,i{k,v\\;Ti-To/2)dv\\dk, (14) 

where 

K{oj,~^,v\\) ^Vi{oj - kv\\)V2{uj ~ [k ~ 2ko)c)Vz{uj,k) (15) 

Vi{uj) ^ iTre-^'t'^) - Rw{^{uj))] 

V2{oj)^e-'^'(^)- ^ (16) 
1/3 (^,fc) = e(.^.2/cofc)g(-fc) 



'^0 



The selective properties of K{u!, k ,v\\) in the longitudinal component of the ve- 
locity I'll are determined by the first cofactor Vi in (^5|). A maximum K{u!, k ,v\\) 
in ti|| at a fixed to is determined by the condition Vi{lu — kv^ g) = max which, as 
a consequence of the properties of the exponential and the Rw functions, cor- 
responds to the frequency Doppler shift condition in the case of the scattering 
from a single particle: 

-ll.o-^- (18) 

The width Aw|| of a maximum K{uj, /c , W||) in wh (that determines the region of 
velocities making the main contribution to the scattered signal at fixed and 
k ) can be estimated from the condition $(a; — kv^ Q ± Aw||) = 1 to be 

At)|| = %/2Aw/fc. (19) 

The selective properties of K(uj, fc in wave numbers k are determined by 
the second cofactor V2 in (p^). A maximum K{uj, k ,v\\) in fc, at a fixed lo, 
is determined by the condition V2(w — (k — 2fco)c) = max which corresponds 
to the condition (analogical to the Volf-Bragg condition for scattering from 
nonstationary spatial harmonic): 

fc = 2fco + -. (20) 
c 

The width Ak of a maximum K{ll!, k ,v\\) in k (that determines the region of 
wave numbers making the main contribution to the scattered signal at a fixed 
uj) can be estimated from the condition ^{oj — (fc — 2fco)c ± Afcc) = 1 to be 

Ak = V2Aw/c. (21) 

The function V3 determines the selective properties oi K(uj, k ,v\\) in the di- 
rection of the wave vectors k and the maximum possible width of the received 
signal spectrum determined by the kinetic function ^(cj, 2fcofc) and antenna fac- 
tor g{—k). 
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5 Discussion 



From ( jT^ ) it follows that the fine structure of the scattered signal spectrum u{uj) 
is determined only by the properties of the function Fn (k ,v\\;T) ( pO| ) and can 
be related to the INPD fii{k ,lt;T) only within the framework of additional 
assumptions about the structure of k ,lf;T). Formally, to determine the 
function fii{k , if; T) require measuring the scattered signal for different wave 
numbers ko of the sounder signal simultaneously. 

To carry out a qualitative comparison with experimental spectra of incoherently- 
scattered signals we use the spectral processing window that repeats the sounder 
signal shape. Their spectra are approximated by Gaussian-like spectra with a 
width equal to their actual spectrum width. According to (p^, the spectrum of 
the received signal is defined by the unknown function Fii{ k ,v\\; T)(p^) (con- 
voluted with a kernel K{ljj, k , For comparison with experimental data, we 
give following simple model of the function Fii( k , T). 

Simple model Assume that /ii(fc,17';T) involves only one peak in the 
range of longitudinal wave vectors and velocities of our interest: 

fii{~k-, If; T) = 5{t - kl)S{lf - vt). 

Fii{ k ,v\\;T) = 5{k — ki)6{k vl — kv\\) . 

This model corresponds to the fact that the medium involves an isolated spatial 
harmonic fii{k ,lf;T) at the wave vector fci and with longitudinal velocity 
vi , the amplitude of which is significantly larger than the amplitudes of spatial 
harmonics close to it. The spectrum of the received signal will then involve also 
only one peak: 

u{co) ^Vi{uj- kiVi)V2{uJ - (fci - 2ko)c)V3{u, fci), (22) 

and the form and width of this peak will be defined by the product Vi{uj — 
kiVi)V2{L0 — (fci — 2fco)c). From the position of the peak in the spectrum 
u(u}), one can determine the unknown wave number fci of the spatial harmonic 
/ii( fc , if; T) that make the main contribution to the scattered signal (|20| ) and 
the longitudinal velocity vi (^8|) corresponding to the moment T . 

This model gives the signal spectrum with only one peak. According to the 
model, the observed presence of several peaks in the real spectrum corresponds 
to existance not only one peculiarity described by the model but a numer of 
them. 

Let us show that after averaging this model gives as the well known expres- 
sion for the IS average spectral power. According traditional approaches lets 
suppose that < fii(k,v,T) >= fio{v){when\k\ > 0). Summarizing over the 
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different realizations gives the produces the function: 



\Vi{uj^ 2koVn ,)|' / \V2{u; - (fc - 2ko)c)fdk 

2 



(23) 



const 



JV3{uj,k)dk J f,o{v)\Vi(iJ~2kov)\'dv 

2 



Ge(i^,2fco)g(-fco) 
8fej;£||(a;,2feo) 



/ /,o(t;)|yi(c^-2fco«)rci^; 



From (g3|) becomes clean, that the average received signal spectral power 
< > is determined as a product of the kinetic function 2fcofc)p 

and maxwellian ions distribution fio{v) convolved with a narrow-band function 
|Vi(a; — 2kov)\'^ (which is defined by the sounding signal and spectral processing 
receiving window spectra) . Qualitatively this solution (^) is close to the average 
spectral power obtainded by the tradditional way (^. 

A numerical simulation have been made for the model containing a number 
of peaks: 



Fii{k, vi\; T) = b{sin{kA + 0(r)) * sin(B{v\\) + 0(T)) - 1), 



(24) 



where A - some constant; B(v) ^ exp{—x^ /v^^)dx] (t){T) - a function with 
uniform distribution of values (we have used random one, having fixed value for 
fixed T ). This expression for Fn has a normal distribution over v\\, uniform 
distribution over fc, and statistical independence of the values for different T (or 
when delay between them exceeds the interval in which (j){T) changes slowly). 
Spectra obtained by substituting this simple model ( p4| ) into obtained equation 
(H), gives us the following single spectrums and their statistical properties, are 
shown at Figure ^. As one can see, single realizations of the spectral power(E- 
G) have the same structure, as an experimental one (Figure]^), the close peak 
width (B) (Figure |[B), and same relation between mean peak appearence (C) 
and average spectral power (D) (Figure ||,C,D). This allows us to suppose that 
the model of the scattered signal single spectrum (^) could be used to describe 
signal properties, and simplified model (|2^) could quilitatively describe the ion 
number phase density behavior. 



6 Conclusion 

In this paper we have suggested an interpretation of separate realizations of in- 
coherently scattered signal spectra. It is based on the radar equation [Berngardt 
and Potekhin,2000] and kinetic theory of ion-acoustic oscillations of a statisti- 
cally homogeneous unmagnetized, coUisionlcss ionospheric plasma with one sort 
of ions [Landafi, 1946]. 

In accordance with the proposed model (^), the main contribution to the 
scattering is made by plasma waves caused by spatial harmonics of ions number 
phase density fii{k ,lf;T), with wave numbers on the order of the double wave 
number of the sounder signal A: w 2fco, for T = Ti — To/2. 



14 




Figure 4: Numerical simulation results for — Ti — 1500K using equations 
(||), ( p^ and simple model (p^. Mean peak amplitude (A), mean peak width(B), 
mean peak appearence(C), average spectral power (D) and three spectral real- 
izations (E-G). 
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It has been shown that the form of the received signal spectrum u{uf) is 
related to the INPD fii{ k ,lf;T). At each frequency uj the value of u{lu) is 
determined by Radon's-like integral Fn ( A: , U| | ; T) on fii{ k ,lf]T) between the 

limits of the velocity components across the wave vector k (^0|). The region 
of wave vectors and longitudinal velocities making contribution to the received 
signal fc ,ti|| is determined by the weight volume K{uj, fc (||). 

So, by changing the transmitting pulse start moment and the moment of its 
receiving one could measure the value of Fn { k ,v\\;T), as function of the time T. 
This allows the INPD fn^k , if; T) diagnostics for different moments including 
delays smaller than theirs lifetime and without statistical averaging of receiving 
signal. Actually, in the case of irregularities lifetime much longer than sounding 
pulses repeating interval T^+i — . , one could measure the Fii( k , W[ [ ; — To/2) 
behavior as function of the time t = Ti — To/2. 

Based on the proposed model in (||) and a Gaussian approximation of the 
spectra of the sounder signal envelope and the receiving window, a qualitative 
comparison of the model with experimental data from the Irkutsk incoherent 
scatter radar was carried out. The comparison showed a quilitative agreement 
for simplified model Jp^) , based on additional assumptions about the properties 
of the function Fn ( fc , V[ [ ; T) . 
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